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. The most distinctive feature of the primordial density inhomogeneities that existed 

\ in the cosmic plasma at the instant of decoupling is their scale invariant power 

spectrum V{k) ~ k over the range k <C H^q of modes with cosmologically large 
CIh' comoving wavelength. We characterize this feature in real space, in terms of their 

2 i correlation function at two points. We show that over cosmologically large comoving 

I distances r ^ H^^ the primordial inhomogeneities were (anti) correlated as /(r) ~ 

— r~^. We revisit the so-called origin of structures problem of the standard cosmology 
at the light of this observation. Our conclusions contradict the current wisdom on 
this issue. 
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Q ■ 1. Tiny temperature anisotropies over the sky in the blackbody spectrum of the cosmic 

Q , microwave background radiation (CMBR) were first detected by the Cosmic Background 



X 
J3 



Explorer (COBE) at the beginning of the 90's [1] and their amphtude and statistical cor- 
relations were precisely measured with angular resolution of a fraction of a degree by the 
Wilkinson Microwave Anisotropy Probe (WMAP) a decade later \^\. Long before, during 
the 70's, these temperature anisotropies had been theoretically predicted |3|, Ij] as the im- 
prints of the primordial seed of tiny density inhomogeneities that must have existed in the 
cosmic plasma at the time of decoupling and from which the present large scale structures 
of the universe would have developed through gravitational collapse during the subsequent 
stage of matter dominated expansion. The measurement with higher precision and better 
spatial resolution of the primordial seed of cosmological density inhomogeneities, as well 
as the analysis and interpretation of its statistical features, is the most crucial enterprise 
in modern cosmology [s], because of the wealth of information that it can disclose on the 
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geometry and content of the universe, as well as on the physical mechanisms that operated 
in it during its early history. 

The available observational data on this random seed of primordial cosmological den- 
sity inhomogeneities is compatible with the assumption of gaussianity (see, for example, 
j^) and, therefore, it is generally accepted that their statistical correlations at two points 
f{x,y) completely describe them. It is also generally assumed, along the fundamental cos- 
mological principle of standard Friedmann-Robertson- Walker (FRW) cosmology, that these 
inhomogeneities are statistically homogeneous and isotropic. Under these assumptions the 
statistical correlations at two points are a function only of the comoving distance between 
them, f{x, y) = f{\x — y\). It is quite common in the literature to represent these correlations 
in terms of their power spectrum in Fourier space (see, for example, [3]), 

r(k) = \ f d'f f{r) e-''-^^ = ^ Hdr f{r) (1) 

The power spectrum has units of volume and is necessarily a non-negative function, 
V{k) > 0. The normalization factor po in this definition is the average cosmic density 
per unit comoving volume. For the sake of simplicity, we use throughout this paper co- 
moving coordinates defined with the scale factor of the universe at decoupling as reference, 

a{teq) = 1. 

2. The most characteristic statistical feature of the primordial seed of cosmological density 
inhomogeneities that existed in the cosmic plasma at the time of decoupling is the scale 
invariance of their power spectrum in the range of wavenumbers with comoving wavelength 
much larger than the comoving size of the FRW causal horizon at that instant H~^^, 

V{k) = 0{k), at A; < Heg. (2) 

The notation 0{kP) used here is borrowed from calculus to denote an infinitesimal function of 
order p. Scale invariance means p = 1, that is, the power spectrum of the primordial density 
inhomogeneities is an infinitesimal of linear order in the range of modes with cosmologically 
large comoving wavelength. Let us remind that the comoving size of the FRW causal horizon 
at decoupling H~g^ is roughly two orders of magnitude shorter than the comoving size Hq^ 
of the causal horizon at present, Hq^ ~ 100 H^^^. 

Scale invariance, which was theoretically predicted as early as the the very existence 
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of the primordial seed of cosmological density inhomogeneities in the radiation dominated 
co.n,.c plasma g ha= been ob.e™tionally co„fi.n,ed by the data collected by the WMAP 

[2j. Nonetheless, this feature ([2]) still poses a major theoretical challenge for the standard 
FRW cosmological framework: widely reported analysis seem to indicate that the power 
spectrum of the primordial density inhomogeneities that could have been generated in the 
radiation dominated cosmic plasma by standard physical mechanisms should be suppressed 
in the cosmological range k <^ H^q much stronger than linearly. Indeed, as V{k) = 
an infinitesimal of, at least, fourth order. 

The theoretical reasoning behind this claim goes as follows. First, it is argued that in the 
framework of standard FRW cosmology causal constrains forbid statistical correlations over 
comoving distances longer than the horizon and, hence, the power spectrum in the range 
k <^ H~g^ can be estimated by Taylor expanding the basis of Fourier functions in ([1]) around 
A; = 0, 



= ^ f dr f{r) (l " ^ + Oiikrf)^ . 
The spatial integral is then moved rightward before the Taylor summation to get. 



(3) 



47r 

Vik) = ^ 



dr /(r) 



Wo 



dr /(r) 



k^ + 0{k^). 



(4) 



As the random density inhomogeneities are constrained to preserve the total mass (energy). 



An / dr r^ fir) = 0, 

JO 

as well as the moment of inertia of the distribution, 



(5) 



An / dr /(r) = 0, 
Jo 



(6) 



it is straightforward to conclude that their power spectrum must be an infinitesimal of, at 
least, fourth order in the cosmological range. 



Vik) = 0{k^) at A; < Heq. (7) 

The discrepancy between the scale invariant power spectrum of the actual primordial cos- 
mological density inhomogeneities ([2]) and the widely reported theoretical prediction ([7]) is 
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commonly known in the literature as the origin of structures problem of standard FRW cos- 
mology. The widely accepted conclusion inferred from this disagreement is that the scale 
invariant power spectrum ([2]) of the primordial cosmological density inhomogeneities is a 
pristine signature of an early period of inflationary cosmology, which cannot have been al- 
tered later on during the subsequent stage of radiation dominated expansion of the universe 
(see, for example, 9|]). In consequence, most efforts in theoretical cosmology during the 
thirty years that have passed since this argument was first laid down have been directed 
towards the implementation of the inflationary paradigm [l^ in a theoretically consistent 
cosmological framework compatible with the available observational data 

3. The motivation of this paper is rather different. We intend to test an implicit and 
unadverted assumption in the reasoning that leads to the statement of the origin of struc- 
tures problem of standard FRW cosmology. In particular, we notice that in order to justify 
expansion (jlj) it is crucial to assume /(r) ~ for r > H~^^, that causal constrains in the 
standard FRW cosmological framework force the two points correlation function to cancel 
or, more precisely, to be strongly suppressed at comoving distances larger than the causal 
horizon. Indeed, if the two points correlation function has compact support it is obviously 
legitimate to perform the spatial integral prior to the Taylor summation. 

The permutation in the order of performing these two operations is not legitimate, 
notwithstanding, if the suppression of the two points correlation function at distances longer 
than the horizon is not strong enough. This can be readily seen from the simple fact that 
the global conservation constrain (I5l), which forces the two points correlation function to be 
suppressed as /(r) ~ o{l/r^) as r ^ oo, guarantees that expression ([T]) is properly defined, 
but it does not guarantee that expression (jl]) is also proper. 

This observation is pertinent because causal constrains in the standard FRW cosmology 
do not necessarily force that the statistical correlations over comoving distances larger than 
the causal horizon do vanish. Actually, the global constrain ([5]) can force /(r) ^ 
for r > H~^^, non-zero statistical correlations in the radiation dominated cosmic plasma 
over comoving distances larger than the causal horizon. This can be readily seen from the 
following simple consideration. We notice that the radiation dominated cosmic plasma before 
decoupling supports the propagation of weak shock waves that travel with sound velocity 
Cs = (in natural units). These waves erase all density fluctuations with comoving 

wavelength shorter than the sound horizon and, hence, force positive statistical correlations 
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between comoving sites separated by shorter distances. That is, /(r) > for r < H~^, 
which is obviously not compatible with the demand that statistical correlations do vanish 
over comoving distances longer than the horizon, since the global conservation constrain ([5]) 
must be fulfilled. 

Let us remark at this point that / dr /(r) = ■^=?- lirrik^o 'P{k) = and, there- 
fore, the global constrain ([5]) is necessarily fulfilled in any causally constrained cosmo logical 
framework. 

4. In order to perform a careful estimation of the limiting behavior of the power spectrum 
([T]) in modes k <^ Heq for two points correlation functions that, in general, have non-compact 
support we separate, 

V{k) =Vshort{k)+Vlong{k), (8) 

the contributions of statistical correlations over comoving distances shorter and larger than 
the causal horizon. 



^ ,,, An f^eq , , ^ , , sin(kr) 
Vshortik) = — drr^ fshortir) — 9 

Po Jo kr 



_ 47r /""^ , 9 „ , , sin(kr) , , 

Viongik) = — drr^ fiong{r) . 10 

On Jh-^ ^ kr 



Air f°° ^ 2 J. ^ ^ sin{kr) 
_ 

We have chosen to denote the projection of the two points correlation function /(r) 
upon the sphere of causally connected comoving sites with the label short, fshortir) = 
/(r)l^gp j^-^i], and the projection upon its complementary volume with the label long, 
fiongir) = fir)le[H-\oo)- As we have noticed above, fiong{r) ~ oil/r^) as r oo for 
the constrain ([5]) to be properly defined. 

The contribution (Q of the two points correlations between causally connected comoving 
sites can be safely estimated in the range of interest k <^ ifgg by Taylor expanding the 
Fourier functions around = 0, as we did in (j4]), Vshortik) = ^ /o ""^ dr fshorti'f') — 

/o dr r^ fshort{r) k'^ + 0{k^). We then use the global constrain ([5]) to write the 
first term in the right hand side of the last expression as a contribution from statisti- 
cal correlations over comoving distances longer than the horizon, /q dr r fshortir) = 
— J^-i dr fiongi^), and get that the total power spectrum ([8]) in the range of very small 
wavenumbers k <^ is. 



6 



V{k) 



— ^drr^ f longer) ' ^ 



3!pg 



/ c??^ r fshortir) 
Jo 



(11) 



The object of interest of our revised analysis is the contribution of the two points correlations 
over long distances, 



V^{k) = — / _^dr r'^ fiong{r) 



sin{kr) 
kr 



(12) 



If these correlations are strongly suppressed this contribution is, not surprisingly, an infinites- 
imal of, at least, second order. More precisely, fiongij') = n > 5, as r oo =^ 
V^lk) = 0{k'^), for k -C Hgg. On the other hand, when the two points correlations at long 
distances are only mildly suppressed. 



fu 



-C pI H-^'+^^ r-(^+^), e e 



(13) 



llongy ) — ~ r\) ^^eq 

where C is some dimensionless constant, their contribution to the power spectrum in modes 
with cosmologically large comoving wavelength is an infinitesimal of order lower than two. 



V,{k) = -47r C H;}^+'^ 



dx 



1 



kH- 



X 



2+e 



sinyx) 



X 



0{k 



(14) 



because as — the integral within the squared brackets tends to a finite negative con- 
stant, — oo < A = linix^o fx dx -^^^ sm{x) _ ^ Actually, the correlations over long 
distances ( IT3i) give a scale invariant contribution to the total power spectrum (JTTl) if e ^ 0, 



V{k) ~ -4 TT C H;J^ Ak- 



An 
Wo 



H- 



dr fshortir) 



e + 0{k' 



(15) 



That is, 



/•oo 

/ dr r"^ f{r) = Q and fir) = 0{l/r^), as r ^oo ^ r{k) = 0{k), fork^Heq 
Jo 

This derivation proves, as a side result, that scale invariant density fluctuations (fT3|) are 
necessarily anticorrelated at long distances, C > 0, because the power spectrum f|T5l) must 
be non-negative everywhere, V{k) > 0. Actually, the global constrain ([5]) imposes that 
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C dr f{r) = /o^- dr r' U^r) ~ C pi H'' /^-i dr r'^ = 
= Jo dr r^ fshort{r) - C pi H'^ = 0, and, therefore, 

C = ^ dr Uortir) = /(O). (16) 

That is, the slope of the scale invariant power spectrum (ITSl) is fixed by the average corre- 
lations within the causal sphere /(O) = 3 H^^ Jq dr fshort{r) and the size of the causal 
horizon at decoupling, V{k) ^ — /(O) -f^^*^ A; for A; <^ ifgg- 

5. This revised analysis naturally lead us to wonder if causal mechanisms operating in 
the cosmic plasma during the radiation dominated expansion of the universe could generate 
mildly suppressed /(r) ~ — energy density (anti) correlations over comoving distances 
longer than the causal horizon. A simple toy model shall convince the reader that this is 
indeed feasible. 

We consider a homogeneous and isotropic fluid that fills the infinite Minkowski space 
and study its response to a local density inhomogeneity generated within it. We describe 
the fluid dynamics within the framework of relativistic hydrodynamics in order to carefully 
account for causality constrains. In this framework the flow equations are 

d^T'^^ix) = 0, (17) 

where 

T^^(s) = [e{x) + P{x)] u^'{x)u''{x) - P(x)^^^ (18) 

is the isentropic energy-momentum tensor and g^'^ = diag{+l, —1, —1, —1) is the Minkowski 
metric tensor. The fluid four-velocity u^{x) = {l,v{x,t)) / \Jl — \v{x, t)^ transforms as a 
Lorentz vector and fulfills the normalization equation u^{x)u^{x) = 1. Finally, e{x) and 
P{x) are, respectively, the fluid internal energy and pressure in its proper local frame. They 
are, obviously, invariant under Lorentz transformations. 

We assume that these two thermodynamic quantities are related by an equation of state 
of the kind. 



P(x) = c1 six), (19) 
where < < 1 is the sound velocity in the fluid. We plug this relationship into (fT8|) 
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T'"^(x) = (l + cl) e{x)u^'{x)u''{x) - cle{x)g^ 
and write down the flow equations 



(20) 



u''{x)u''{x)df,e{x) + e{x) [u''{x)d^u''{x) + u''{x)d^u>'{x)] 



[l + cl 



-g^''d,e{x) = 0. (21) 



We find convenient to formulate them as follows. First, we multiply (|2Ti) by Ui,{x) to obtain 
the scalar equation 



u''{x)d^e{x) + (l + cfj e{x)df,u^'{x) = 0, 
which we plug into fl2T]) to find 



(22) 



e{x)u^{x)d^v!'{x) + 



l + cl 



u^'{x)u''{x) - g^""] d^e{x) = 0. 



(23) 



We now specialize these fully covariant equations to the preferred frame located at the origin 
of a spherically symmetric fiow, e{x,t) = e{r,t), v{x,t) = v{r,t)f. As we said, we assume 
that this fiow is the response to a local density inhomogeneity generated at the origin. We 
find 



dte{r, t) + f (r, t)(9,e(r, t) + (l + cf) e(r, t) 



v{r,t)dtv{r,t) + drv{r,t) 2 

^— + -v{r, t) 

1 — t>^(r, t) r 



€{r, t) ' i ^ ' + -r^, [v{r, me{r, t) + 9,e(r, t)] = 0. 



1 — t>^(r, t) 



1 + d 



This system of hyperbolic equations can be numerically solved, in general, by the method of 
the characteristics [l^. Its solutions show [l^ that the radial fiow can extend well beyond 
the causal light-cone of the origin. This can be easily seen, in particular, when the fiow is 
stationary: 



v{r)dre{r) + (l + cfj e{r) 



1 — v^yr) r 



(24) 
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e(r)^P^ + Th^Mr) = 0. (25) 

l — v^[r)l + cj 

It is then straightforward to find 

\ cf y 1 — f^(r) r 

According to this equation the fiow velocity v{r) everywhere beyond the hght-cone is fixed 
by the fiow velocity vr at the light-cone surface. In particular, when f ^ <^ < 1 the fiow 



is roughly incompressible 



vir) ~ vn (-) (27) 
The energy density associated to this fiow is given, in this asymptotic limit, by 

T°°(r) ^ e{r) ^ sn e'^^'-'^) ^ ■ [l - ^4^) • (28) 

That is, the response of a fiuid to a local density inhomogeneity can generate scale invariant 
density (anti) correlations over scales well beyond its light-cone. 

These conclusions challenge the prevalent wisdom in cosmology, according to which the 
scale invariant power spectrum of the primordial cosmological density inhomogeneities is a 
pristine signature of the very early history of the universe that could have not been altered 
during the subsequent FRW standard expansion. 
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